op 



Cell Research (2013) 23:1106-1118. gm^ 
©2013 IBCB.SIBS, CAS All rights reserved 1001-0602/13 HiP 



ORIGINAL ARTICLE www.nature.com/cr 



The gating charge pathway of an epilepsy-associated 
potassium channel accommodates chemical ligands 
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Voltage-gated potassium (Kv) channels derive their voltage sensitivity from movement of gating charges in volt- 
age-sensor domains (VSDs). The gating charges translocate through a physical pathway in the VSD to open or close 
the channel. Previous studies showed that the gating charge pathways of Shaker and Kvl. 2-2.1 chimeric channels are 
occluded, forming the structural basis for the focused electric field and gating charge transfer center. Here, we show 
that the gating charge pathway of the voltage-gated KCNQ2 potassium channel, activity reduction of which causes 
epilepsy, can accommodate various small molecule ligands. Combining mutagenesis, molecular simulation and elec- 
trophysiological recording, a binding model for the probe activator, ztz240, in the gating charge pathway was defined. 
This information was used to estabUsh a docking-based virtual screening assay targeting the defined Ugand-binding 
pocket. Nine activators with five new chemotypes were identified, and in vivo experiments showed that three Ugands 
binding to the gating charge pathway exhibit significant anti-epilepsy activity. Identification of various novel acti- 
vators by virtual screening targeting the pocket supports the presence of a ligand-binding site in the gating charge 
pathway. The capability of the gating charge pathway to accommodate small molecule Ugands offers new insights 
into the gating charge pathway of the therapeutically relevant KCNQ2 channel. 
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Introduction 

Voltage-gated potassium (Kv) channels represent the 
largest group of potassium channels and play crucial 
roles in a wide range of physiological processes [1, 2]. 
Kv channels constitute potential drug targets for the treat- 
ment of diverse diseases ranging from cancers to meta- 
bolic, cardiovascular and neurological disorders [3]. The 
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six transmembrane segments of Kv channels are grouped 
into the voltage-sensor domain (VSD) (S1-S4) and the 
pore domain (S5-P-S6). The voltage-sensing function of 
the VSD is conferred by charged residues, mostly argi- 
nines in S4. To open or close the channel, these charged 
residues translocate between the internal and external 
water-accessible crevices in the VSD through a physi- 
cal gating charge pathway in response to membrane 
potential changes [4]. Studies on Shaker and Kvl. 2-2.1 
chimeric channels suggest that the pathway is occluded 
in the middle of the VSD [5-1 1]. A conserved phenylala- 
nine in the middle of S2 (F290 of the Shaker channel and 
F233 of the Kvl. 2-2.1 chimera channel), together with 
the nearby arginines of S4, was proposed to form the oc- 
cluded site [5, 7, 12]. The phenylalanine separates the 
extracellular and intracellular hydrated crevices of the 
VSD and forms the charge-transfer center that catalyzes 
movement of the gating charges [5, 13]. Substitutions of 
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certain arginines of Shaker S4 by smaller amino acids 
create an omega pore in the mutated VSD, which is per- 
meable to some cations (e.g., K^, Li^ and Cs^) [7, 12, 14]. 
The omega pore has been demonstrated to be the same 
as the gating charge pathway [5], which has the proper 
size to accommodate an arginine side chain [15]. The 
occluded pathway is consistent with the focused electric 
field theory, i.e., the transmembrane electric field is fo- 
cused on the occluded site by the internal and external 
water crevices [7, 16]. In addition, several previous stud- 
ies showed that the arginines along S4 are accessible to 
methanethiosulfonate reagents [8, 17, 18], implying that 
small molecules might bind to this site. The compounds 
ICA-27243 and NH29 were also suggested to act on the 
VSD [19, 20]. These studies shed light on the central 
structures of the gating charge pathways of the Shaker 
channel and the Kv 1.2-2.1 chimeric channel, and raised 
the possibility of the gating charge pathway as a thera- 
peutic target [5, 13]. 

KCNQ2 potassium channel (Kv7.2) is a low-threshold 
Kv channel. Opening of KCNQ2 channel or heteromul- 
timeric KCNQ2 complexes below threshold will inhibit 
initiation of action potentials, and thus KCNQ2 plays a 
fundamental role in controlling neuronal excitability [21]. 
Mutations of human KCNQ2 genes that result in reduc- 
tion or loss of channel activity cause benign familial 
neonatal convulsions. Activators capable of augmenting 
KCNQ2 function have been proven effective in treatment 
of human epilepsy [22]. In this study, we have identified 
an activator-binding pocket in the gating charge path- 
way of KCNQ2 by comprehensively employing homol- 
ogy modeling, molecular docking, molecular dynamics 
(MD) simulation, mutagenesis and electrophysiological 
determination. The small molecular activator that we 
recently discovered, ztz240 [23], was used as a probe for 
determining the binding model of ligands in the gating 
charge pathway. zt240 and two of the nine newly identi- 
fied activators targeting the gating charge pathway ex- 
hibited comparable anticonvulsant activity as retigabine, 
a KCNQ activator approved as an anti-epilepsy drug 
in 2011 [24]. Our results offer insights into the gating 
charge pathway of KCNQ2 channel and demonstrate that 
the pathway can serve as a drug target for development 
of new therapeutics. 

Results 

Essential roles of the conserved phenylalanine (F137) 
for ztz240 activity 

Compound ztz240 is a KCNQ2 activator identified 
recently [23]. Major effects of ztz240 on KCNQ2 chan- 
nel include increasing outward current amplitude, left- 
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shifting the voltage-dependent activation curve {G-V 
curve) and dramatic slowing of deactivation (Figure 1). 
The phenylalanine forming the occluded site in Shaker 
and Kv 1.2-2.1 chimeric channels is highly conserved in 
the family of Kv channels with the degree of conserva- 
tion of ~99% [25]. The corresponding conserved phe- 
nylalanine in KCNQ2 is F137, which is located around 
the suggested occluded site of gating charge pathway 
(Figure lA and IB; Supplementary information. Figure 
SI A). The activator ztz240 is mainly composed of two 
aromatic rings with hydrophobic properties linked via an 
amide bridge; we therefore speculated that it may form 
hydrophobic interactions with F137 if the compound 
could access this residue. To test this idea, we performed 
mutagenesis on the conserved residue. Indeed, muta- 
tion of KCNQ2 F 137 to alanine (F137A) dramatically 
reduced ztz240 activity (Figure IC-IF). The right-shifted 
dose-response curve of KCNQ2 F137A fiirther supports 
the critical role of F 137 for the activity of ztz240 (Fig- 
ure IG). For the wild-type channel, the dose-response 
curve was a well-fit sigmoidal-shaped curve with an EC50 
value of 2.81 ^iM. In contrast, the dose-response curve of 
KCNQ2 F137A could not be fit well, as potentiation is 
rather low even at higher concentrations. 

One possible cause of the reduced sensitivity is that 
the transition into or out of open state of KCNQ2 chan- 
nel is affected by F137A. However, negligible changing 
of biophysical parameters argues against F137A perturb- 
ing "state transition". The tested biophysical parameters 
of KCNQ2 F137A, such as current density, half maxi- 
mal activation voltage (Fiq), time constants of activa- 
tion (closed to open transition) and deactivation (open 
to closed transition), are very similar to those of wild- 
type channel (Supplementary information. Figure S2A). 
In addition, F137A did not reduce the activities of other 
KCNQ2 activators. Retigabine, an activator recognizing 
a tryptophan in the S5 segment [26, 27], was tested first. 
KCNQ2 F137A exhibited similar sensitivity to 10 ^M 
retigabine to wild-type channel (Supplementary informa- 
tion. Figure S2B). Next, ZnPy, another activator potenti- 
ating KCNQ channels through different mechanism from 
retigabine [28], was found to potentiate wild-type and 
F137A KCNQ2 channels in comparable scales (Supple- 
mentary information. Figure S2C). NH29, the third acti- 
vator reported to target the external part of the VSD [20], 
also potentiated wild-type and F137A KCNQ2 channels 
in a similar manner (Supplementary information. Figure 
S2D). These data indicate that F137 is specifically criti- 
cal for activity of ztz240. 

Identification of more critical residues in the KCNQ2 VSD 
F137 is located at the middle of the S2 segment of 
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Figure 1 F137 of KCNQ2 is critical forztz240 activity. (A) The crystal structure of ttie VSD of Kv1 .2-2.1 chimeric channel (PDB 
code: 2R9R). The conserved phenylalanine in S2 (F233) is represented by spheres. Charged residues in S4 are displayed 
as sticks. The blue shadow indicates the occluded site. (B) Sequence alignment for S2. The conserved residues are high- 
lighted in red. The conserved phenylalanine in S2 is indicated by an arrow. (C) Representative traces of wild-type and F137A 
KCNQ2 channels with or without 10 nM ztz240. The inset shows the chemical structure of ztz240. (D) 10 nM ztz240 potenti- 
ated outward current amplitude of wild-type, but not F137A KCNQ2 channels. The measured currents were elicited by +50 mV. 
(E) G-V curves of wild-type and F137A KCNQ2 channels with or without ztz240. For wild-type KCNQ2, 10 \iU ztz240 left- 
shifted the Vi,2 dramatically, whereas the change for KCNQ2 F137A was negligible. (F) Histograms show the time constants 
of deactivation of wild-type and F137A KCNQ2 channels with or without ztz240. Deactivation phases that follow +50 mV 
depolarization were analyzed. For wild type, the deactivation phase was fit by a single exponential function before application 
of ztz240. However, the prolonged deactivation phase after application of ztz240 was fit by a biexponential function with fast 
(TaUfas,) and slow (TaUs^^) time constants. For KCNQ2 F137A, both the deactivation phases with and without ztz240 were fit 
well by a single exponential function. (G) Dose-response curve of ztz240 on wild-type and F137A KCNQ2 channels. The test 
potential is +50 mV. 



KCNQ2 VSD. We designed more mutants for the resi- 
dues lining the VSD. Among the functional mutations, 
E130A, I134A, G138A, R207W and R2 1 OA dramatically 
decreased the potentiation activity of ztz240 in all three 



major aspects. They reduced the increase of outward cur- 
rent amplitude ratio (///q), and prevented or attenuated 
the left-shifting of the G-F curve (AFi^) and the slowing 
of deactivation (Figure 2). The almost overlapping tail 
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currents of these mutants with and without ztz240 clearly 
show the decreased effects on deactivation (Supplemen- 
tary information, Table SI). In contrast, mutations S105A 
and M174A attenuated the left-shifting of G-F curve but 
did not decrease ///(, ratio or prevent the slowing of deac- 
tivation (Figure 2 and Supplementary information. Table 
SI). Identification of multiple essential residues fiarther 
indicates the importance of the VSD for ztz240 activity. 

Structural model of ztz240 binding to the KCNQ2 VSD 

As the experimental structures of KCNQ2 and its 
complexes with other ligands have not been determined, 
it is challenging to construct a structural model for the 
ztz240-KCNQ2 interaction. Accordingly, a hierarchical 
strategy was used to construct the structural model by 
comprehensively employing homology modeling, mo- 
lecular docking and MD simulation in conjunction with 
mutagenesis and electrophysiological determinations. 
The mutagenesis result suggests that the potential bind- 
ing pocket of ztz240 is mostly located in the VSD rather 
than in other domains (Figure 1). Further electrophysi- 
ological experiments support that ztz240 binds to the 
open-state VSD (Supplementary information. Figure S3). 
Thus, we first built up a three-dimensional (3D) model 
for the transmembrane domain of KCNQ2 based on the 
structural information of the open-state Kvl.2 channel 
[29-31] by using Discovery studio 2.6 (Supplementary 
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information. Figure SIB). 

Second, we constructed the binding model for ztz240 
to KCNQ2. We defined a potential binding pocket for 
ztz240 according to the mutagenesis and electrophysi- 
ological data, which includes F137 and other residues 
essential for ztz240 to potentiate the channel (Figure 2 
and Supplementary information. Table SI). Homology 
modeling could not predict proper conformations for the 
residues interacting with a ligand. Therefore, Induced Fit, 
a flexible docking program that considers the flexibilities 
of both ligand and side chains of a target protein, was 
applied to simulate the binding of ztz240 to KCNQ2. 
However, the current docking methods could not deter- 
mine the orientation of ztz240 entering into the bind- 
ing pocket, indicating that either end may enter into the 
pocket first (Supplementary information. Figure S4A). 
The orientation of ztz240 in the VSD was determined 
with the assistance of two derivatives of the compound 
(named YGPl and YGP2) that have the fluorine and the 
chlorine substituted by hydrogen atoms, respectively. 
Electrophysiological recording showed that YGPl still 
exhibited potentiation activity on KCNQ2 as ztz240 did; 
however, the potentiation activity of YGP2 was dramati- 
cally reduced (Supplementary information. Figure S4B 
and S4C). This result indicates that ztz240 may bind to 
the VSD via the configuration with chloropyridine form- 
ing the bottom of the pocket and the fluorophenyl located 
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Figure 2 Identification of more critical residues in KCNQ2 VSD for ztz240 activity. (A) Representative traces of identified criti- 
cal mutants with and without 10 nM ztz240. (B) Effects of 10 )jM ztz240 on outward current amplitude of VSD mutants (n > 
3). Each mutation site is indicated based on the predicted transmembrane regions. The dashed line indicates a potentiation 
level of 10 |xM ztz240 on the wild-type KCNQ2 channel. The test potential is +50 mV. (C) Effects of 10 |xM ztz240 on Vy2 of 
VSD mutants [n > 3). AVy2 indicates the change of after application of ztz240. AVy2 = Vy2 in the presence of ztz240 - Vy2 
in control. 
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near the pocket entrance (Supplementary information, 
Figure S4B). 

Third, the ztz240-KCNQ2 binding model derived 
from the docking simulation was further optimized by 
MD simulation. The MD simulation may overcome the 
shortage of docking simulation that could not consider 
the flexibility of the protein backbone [32]. The model 
of the tetramer of ztz240/KCNQ2 complex was embed- 
ded into a patch of l-palmitoyl-2-oleoyl-sn-glycero-3- 
phosphocholine (POPC) bilayer and then subjected to 
a 200-ns MD simulation. The ligand-VSD interactions 
suggested by the optimized binding model are shown in 
Figure 3. ztz240 hydrophobically interacts with 1134 and 
F137; the nitrogen atom of the linker forms hydrogen 
bonds with the side chain of E130; the pyridine ring is 
involved in a CH-;i interaction with the a-hydrogen of 
R207; the chlorine atom forms an electrostatic interac- 
tion with the guanidyl of R210 and a halogen bond with 
the carbonyl oxygen of R207, respectively; the pyridine 
nitrogen forms another electrostatic interaction with the 
guanidyl of R2 1 0. Residues that have been identified as 
critical for the potentiation activity of ztz240 in the mu- 
tagenesis studies are directly involved in the ztz240-VSD 
interactions, which explains well the mutagenesis data 
(Figure 2 and Supplementary information. Table SI). 

Most of the crucial residues lining the binding pocket 
are located in S2 and S4, such as E130, 1134, F137, 
G138, R207 and R210, while very few residues in SI 
and S3 affect the compound, such as SI 05 and Ml 74 



(Figure 3), which are not conserved in KCNQ isoforms. 
To further evaluate the roles of residues in the external 
segments of S 1 and S3 in ztz240 activation, the remain- 
ing residues except those listed in Figure 2B and 2C 
were examined by alanine scanning, and none of them 
significantly affected ztz240 activity (Supplementary 
information. Table SI). These additional mutations and 
electrophysiological determinations retrospectively vali- 
date again our ztz240-KCNQ2 binding model. 

The model of the ztz240-KCNQ2 interaction clearly 
describes a broad pocket in the VSD, extending from the 
extracellular mouth of the VSD to R210 (Figure 3C). 
This possible ztz240-binding pocket in KCNQ2 partially 
overlaps with the corresponding regions of the gating 
charge pathways of Shaker and Kvl. 2-2.1 chimeric chan- 
nels, though it is more expansive and deeper (Figures 
lA and 3C) [5-12, 14, 15]. This result is consistent with 
previous studies indicating that the VSD of Kv channels 
might be draggable sites [19, 20]. Our study further re- 
vealed that the gating charge pathway of KCNQ2 might 
have a new functionality, serving as a direct binding 
pocket for existing and newly designed activators such as 
ztz240 and new compounds discovered in this study (see 
below). 

Discovering new activators targeting the ztz240-binding 
pocket 

To discover new KCNQ2 activators and to fiirther ver- 
ify the possible activator-binding pocket, we carried out 
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Figure 3 The binding model of ztz240 with VSD. (A, B) General (A) and detailed (B) views of the interactions between ztz240 
and VSD. ztz240 and the residues are represented as ball-and-sticks. Carbon, oxygen, nitrogen, chlorine, fluorine and hydro- 
gen atoms in ztz240 are colored in light blue, red, dark blue, green, orange and white, respectively. For clarity, only few key 
hydrogen atoms are shown. Hydrogen bonds or electrostatic interactions are indicated by red dashes. The CH-k interaction 
and halogen bond are indicated by black and green dashes, respectively. (C) The structure of the binding pocket extracted 
from the structural models of ligand-VSD complexes. For clarity, the ligands were removed. The surface of the pocket is 
shown in yellow. Three cross-sections of the pocket from the external view are shown on the right side. Residues are shown 
as sticks. 'Out' and 'In' indicate the external and the internal sides of cell membrane, respectively. 
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a structure-based virtual screen targeting the identified 
pocket, guided by the notion that if the identified bind- 
ing pocket is reasonable, novel activators recognizing the 
identified ztz240-binding pocket could be discovered. 
The SPECS database containing ~200 000 chemicals was 
screened by docking approach. Twenty-five hits (YGOO 1 - 
YG025) were selected for bioassay. Nine compounds be- 
longing to five different chemotypes displayed significant 
potentiation effects on the outward current of KCNQ2 
channel (Figure 4). The potentiation ranges from 1.28 ± 
0.06 to 3.56 ± 0.23 at 30 |iM (Figure 4A). EC50 values of 
the three most potent active compounds, YG002, YG007 
and YG025, were further revealed to be 4.64 ± 0.14 ^M, 
1.60 ± 0.16 |iM and 5.92 ± 0.21 |iM, respectively (Fig- 
ure 4B). The i^M-level EC50 is similar to the approved 
drug retigabine [27]. The scaffolds of the nine identified 
activators can be described as two aryl groups connected 
by a functional linker The linkers can be divided into 



five chemotypes, thiourea (YG002 and YG014), hy- 
drazide (YG006, YG009 and YG025), urea (YG007 and 
YG018), amide (YG021) and amine (YG023). YG021 
is an analogue of the reported activator ICA-27243 [33], 
but the other four chemotypes have never been reported 
before as KCNQ activators. The key residues involved in 
the interaction of ztz240 with the VSD, such as F 137 and 
R210, are still essential for the newly identified activa- 
tors (Figure 4D). Consistently, both F137A and R210A 
mutations were found to differentially affect the activity 
of most of these newly identified ligands (Supplementary 
information. Figure S5). This result demonstrates again 
that the gating charge pathway of KCNQ2 may accom- 
modate various activators. 

Anti-epilepsy activity of the activators targeting the gat- 
ing charge pathway 

A primary anti-epilepsy screening assay was per- 




YGO2I YG023 YG025 



Figure 4 Structure-based virtual screening identifies activators with diverse chemotypes. (A) The increase of the outward 
current amplitude in the presence of compound as indicated. The test potential is -10 mV. The dashed line indicates a poten- 
tiation level of 1 (i.e., no potentiation). (B) Dose-response curves of YG002, YG025 and YG027 (n > 3 for each data point). (C) 
Chemical structures of the discovered activators. The fragments oriented toward the intracellular end of the VSD are high- 
lighted in orange. (D) General view of superimposed binding models of nine identified activators with the VSD. 
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Table 1 Effects of ztz240, YG007 and YG018 on mouse MES- 
induced seizure model 



Group Administration Episodes of seizure Protection (%) 



YG007 


40 mg/kg 


1/10*" 


90*" 


YG018 


40 mg/kg 


1/10*" 


90*** 


ztz240 


20 mg/kg 


0/8*** 


100*" 


Retigabine 


20 mg/kg 


0/8*" 


100*" 


Vehicle 


0.2 ml/10 g 


10/10 


0 



Abbreviation: MES, maximal electroshock. P < 0.001, significantly 
different from vehicle, Fisher's exact probability test. 



formed to evaluate the activity of these compounds in 
mice using the maximal electroshock (MES)-induced 
seizure model. Among the tested activators, ztz240, 
YG007 and YG018 exhibited > 50% protection rate and 
thus were selected for further investigation in MES- 
induced and pentylenetetrazol (PTZ)-induced seizure 
models (Supplementary infonnation, Table S2). Retiga- 
bine was set as the positive control. All three compounds 
exhibited excellent anti-epilepsy activity in both models. 
In the MES-induced seizure model, a single dose admin- 
istration of the three compounds significantly prevented 
the tonic hind limb extension. The protection rates (100%, 
90% and 90%) are comparable with that of retigabine 
(Table 1). In the PTZ-induced seizure model, the three 
compounds exhibited comparable anti-epilepsy activity 
as retigabine. A single dose of these compounds not only 
significantly delayed the clonus latency but also reduced 
the incidence of generalized tonic-clonic seizure (GTCS) 
and mortality (Table 2). 

Discussion 

For many drug targets, such as GPCRs, nuclear recep- 
tors and ligand-gated ion channels, their functions are 
modulated by endogenous ligands. Accordingly, many 
drugs bind receptors at sites where the natural molecular 
activators (i.e., endogenous ligands) act [34]. In contrast, 
the binding site of Kv channel activators is more elusive. 



as the activation of Kv channels in the physiological 
environment is modulated by membrane potential rather 
than endogenous ligands. Elucidation of the binding of 
activators with Kv channels is essential both for under- 
standing the mechanism of channel activation and for the 
development of therapeutics. To date, limited informa- 
tion about the binding site of activators of Kv channels 
is available. The present study provides evidence for the 
presence of a ligand-binding site in the gating charge 
pathway of KCNQ2 channel. The identified activator- 
binding site lying deeply inside the VSD is largely differ- 
ent from postulated sites reported previously, which are 
either located at the pore region or a superficial part of 
the VSD [20, 26, 27]. The newly identified site encases 
the activators through multiple transmembrane segments. 

To characterize a ligand-binding site, crystallography 
of the ligand/protein complex is the most direct method. 
However, membrane proteins are notoriously difficult to 
purify and crystallize. Very few Kv channels have been 
crystallized successfully [35-37], and no drug/Kv com- 
plex structure has been determined. Radioactive ligand- 
binding assays or photo-labeling experiments have not 
answered the question either. As a primary strategy for 
exploring the action sites of activators, mutagenesis 
experiments have been effective in identifying the key 
molecular determinants for potentiation on hERG and 
KCNQ channels [38-42]. However, the limitation of 
mutagenesis is the inability to exclude a binding site 
or a key site that is only essential for the activity of the 
compound, but not for binding. Our study provides an 
alternative way to distinguish binding sites from other 
key sites for channel modulators. We speculated that new 
activators would be found based on the 'binding site' 
by structure-based drug discovery if the 'binding site' is 
true. Structure-based virtual screening has been widely 
applied in drug discovery. A virtual screen of ~200 000 
compounds in the SPECS database targeting this bind- 
ing site was performed. Twenty-five hits were selected 
for biological validation and nine showed potentiation 
on KCNQ2 channel in the voltage clamp assay, which 
is considered as the 'gold standard' for the study of ion 



Table 2 Effects of ztz240, YG007 and YG018 on mouse PTZ-induced seizure model 


Group 


Administration 


Clonus latency (s) 


Incidence of GTCS 


Mortality 


YG007 


40 mg/kg 


161 ±42* 


3/10* 


5/10* 


YG018 


40 mg/kg 


191 ±27*" 


1/10" 


2/10"* 


ztz240 


20 mg/kg 


147 ±35* 


0/10*" 


o/io"* 


Retigabine 


20 mg/kg 


181 ±40"* 


0/10*" 


2/10"* 


Vehicle 


0.2 ml/10 g 


116±26 


8/10 


10/10 



Abbreviations: GTCS, generalized tonic-clonic seizure; PTZ, pentylenetetrazol. P < 0.05, P < 0.01, P < 0.001, significantly different from ve- 
hicle. Fisher's exact probability test. 



Cell Research | Vol 23 No 9 1 September 2013 



Ping Lief a/. @ 



channels. The hit rate of our virtual screening is 36% 
(9/25), which is much higher than the hit rate (0. 1 %) of 
a cell-based high throughput screening for discovering 
KCNQ2 activators (http://pubchem.ncbi.nlm.nih.gov). In 
addition, the fact that all identified novel activators rec- 
ognized the same site further supports the presence of the 
ligand-binding site. It should be noted that the complex 
model was built based on a homology model of KCNQ2 
channel, which may contain several inherent limitations. 
There is no representation of cation-;: interaction and 
halogen bond in empirical force fields. Molecular dock- 
ing and MD simulation applied in the current work could 
not fully characterize these interactions. These particular 
interactions were proposed based on geometric consid- 
erations. Further crystallography or combined quantum 
mechanics/molecular mechanics (QM/MM) studies are 
necessary to fiilly elucidate the binding mode of ztz240 
in the binding pocket. 

The occluded gating charge pathways of Shaker and 
Kvl. 2-2.1 chimeric channels focus the electric field on 
a narrow region and form the charge-transfer center, 
both of which were suggested to be important for the 
mechanism of voltage gating of Kv channels. This study 
shows a new functionality of the gating charge pathway 
of KCNQ2, which serves as a direct binding pocket for 
activators. The broad pocket in the VSD extends from 
the extracellular mouth of the VSD to R210 (Figure 3C) 
and allows chemical ligands to enter and reside in the 
space beneath R207. The volume of this space is ~170 
A^, which is large enough to accommodate the aromatic 
rings of ztz240 and our newly identified ligands (moieties 
highlighted in orange in Figure 4D) with sizes of 90 to 
140 A\ The ability of the KCNQ2 VSD to bind chemi- 
cal ligands suggests that the gating charge pathway of 
the channel might be structurally different from other Kv 
channels, such as the Kvl. 2-2.1 chimera. The phenylala- 
nine (F233) is the key residue of the charge-transfer cen- 
ter in the Kvl. 2-2.1 chimera. When it was substituted by 
19 other natural amino acids, the voltage dependence of 
most mutants was largely affected [13]. However, analo- 
gous mutations F137X (X represents the other 19 natural 
amino acids) cause much less influence on the voltage 
dependence of KCNQ2 channel gating (Supplementary 
information. Figure S6). First, all 19 KCNQ2 mutants 
are functional, whereas only 15 Kvl. 2-2.1 mutants are 
functional. When the phenylalanine is mutated to aspartic 
acid (D), lysine (K) or arginine (R), the KCNQ2 chan- 
nels are still fiinctional while the Kvl. 2-2.1 channels lose 
functionality, suggesting that the phenylalanine (F137) 
in KCNQ2 is more tolerant to changes. Second, effects 
of the mutations on G-F curves are different. Among 
the 15 functional Kvl. 2-2.1 mutants, 13 exhibit largely 



right-shifted G-F curves with positive and the right- 
shifting effects are > 100 mV for some mutants, such as 
F233A [13]. Only two mutations, tryptophan (W) and 
tyrosine (Y), exhibit currents near wild-type levels with 
negative Vyj, which is consistent with the notion that a 
rigid cyclic side chain is important at this position [13]. 
In contrast, for KCNQ2, only six mutants display posi- 
tive values and the right-shifting is < 100 mV for all 
the mutations. Furthermore, besides F137W and F137Y, 
F137A, F137V, F137I and F137L also exhibit similar Vy2 
values to that of the wild type, suggesting that the cyclic 
side chain at F137 position is not as important as in the 
Kvl. 2-2.1 chimera. However, for some specific substi- 
tutions, the mutations lead to similar kinetic effects for 
the two channels. For instance, both Kvl. 2-2.1 F233W 
and KCNQ2 F137W exhibit a slower rate of closure 
than wild-type channels (Supplementary information. 
Figure S6C and S6D) [13]. The discordant phenotypes 
of the F137X and F233X mutants support the idea that 
gating charge pathways are different at this particular 
position. Other explanations for discordant phenotypes 
between KCNQ2 and Kvl. 2-2.1 chimeric channels are 
possible. Interestingly, the present concept of diverse 
gating charge pathways is in agreement with the observa- 
tion that the voltage-gated proton channel Hvl forms an 
internal water wire penetrating the whole VSD [43-45]. 
The Hvl channel only contains four transmembrane seg- 
ments homologous to the VSD of Kv channels, but lacks 
the classic pore domain of other voltage-gated ion chan- 
nels. 

Although there is no available atomic-resolution struc- 
ture of a Kv channel in the resting state, the first four 
arginines (R1-R4) of Kvl. 2 are generally thought to stay 
in a 'down' position in the resting state compared with 
the open state [46, 47]. A recent long-time all-atom MD 
simulation also supports the downward movement of 
arginines [5]. The aromatic moieties of ztz240 and other 
identified activators reside in the gating charge pathway 
by filling the space between R207 and R210, thereby 
preventing the arginine from moving down to its resting 
state (Figure 5 and Supplementary information, Movie 
SI). The capability of the KCNQ2 VSD to accommodate 
chemical ligands allows the possibility to upregulate 
Kv channel function and treat diseases through acting 
at the gating charge pathway directly. This notion was 
examined by in vivo experiments. Anti-epilepsy effects 
of ztz240 and the nine newly identified activators were 
tested in the mouse MES-induced seizure model. Among 
them, ztz240, YG007 and YG018 exhibit anti-epilepsy 
activity and protection rates comparable to that of reti- 
gabine (Table 1). Different from ztz240, retigabine acti- 
vates KCNQ2-5 through a crucial tryptophan residue in 
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Figure 5 The gating cliarge pathway of KCNQ2 channel accommodates chemical ligands. (A) Details of how activators inter- 
act with and affect R207 and R210 in open state, using ztz240 and YG002 as examples. The activators fill the space between 
R207 and R210. The dashed lines indicate the interactions between the activators and R210. (B) The cartoon shows how 
a chemical ligand acts on the gating charge pathway by residing in the gating charge pathway and hindering the arginine to 
move to its down state. 



S5 [26, 27]. As KCNQ channels are accessible to various 
types and combinations of pharmacological regulation 
[28], ztz240 and retigabine may lead to a synergistic ef- 
fect, which might be beneficial for treatment of epilepsy. 
In summary, our study provides new insights into the 
structures and functions of gating charge pathways of 
Kv channels. The discovery of new modulators targeting 
the gating charge pathway by virtual screening may have 
general implications for other voltage-gated ion chan- 
nels. 

Materials and Methods 

Mutagenesis 

Starting from KCNQ cDNAs, the KCNQ2 point mutants were 
constructed by recombinant PCR and verified by sequencing. 

Cell culture and transient transfection 

Chinese Hamster Ovary (CHO) cells were grown in 50/50 
DMEM/F12 (Cellgro, Manassas, VA, USA) with 10% fetal bovine 
serum (FBS) and 2 mM L-glutamine (Gibco, Carlsbad, CA, USA). 
To express ion channels, cells were split at 24 h before transfec- 
tion, plated in 60-mm dishes and transfected with Lipofectamine 
2000™ reagent (Invitrogen, Carlsbad, CA, USA), according to 
the manufacturer's instruction. At 24 h after transfection, cells 
were spht and replated onto coverslips coated with poly-L-lysine 
(Sigma-Aldrich, St Louis, MO, USA). A plasmid for cDNA of 
GFP (Amaxa, Gaithersburg, MD, USA) was co-transfected to aid 
the identification of transfected cells by fluorescence microscopy. 



Electrophysiological recording 

Whole-cell voltage clamp recording was carried out at room 
temperature in CHO cells with an Axopatch-200B amplifier 
(Molecular Devices, Sunnyvale, CA, USA). The electrodes were 
pulled from borosilicate glass capillaries (World Precision Instru- 
ments, Sarasota, FL, USA). When filled with the intracellular so- 
lution, the electrodes had resistances of 3-5 MQ.. Pipette solution 
contained: KCl 145 mM, MgClj 1 mM, EGTA 5 mM, HEPES 10 
mM and MgATP 5 mM (pH 7.3 with KOH). During the recording, 
constant perfusion of extracellular solution was maintained using 
a BPS perfusion system (ALA scientific Instruments, Westburg, 
NY, USA). Extracellular solution contained: NaCl 140 mM, KCl 3 
mM, CaCl, 2 mM, MgCl^ 1.5 mM, HEPES 10 mM and glucose 10 
mM (pH 7.4 with NaOH). Signals were filtered at 1 kHz, and digi- 
tized using a DigiData 1440 with pClamp 9.2 software (Molecular 
Devices). Series resistance was compensated by 60%-80%. 

Mouse MES- and PTZ-induced seizure models 

All animal procedures were performed in accordance with the 
National Institutes of Health Guide for the Care and Use of Labo- 
ratory Animals, under protocols approved and strictly followed by 
the Institutional Animal Care and Use Committees (lACUC). The 
lACUC checked all protocols and approved this study. Male KM 
mice, weighing 1 8-22 g, were obtained from Shanghai SLAC Lab- 
oratory Animal CO Ltd. For the MES-induced seizure model, the 
shock level was set at 160 V and the duration was set at 5.4 s with 
continuous wave under the eighth configuration, using a physi- 
ological and pharmacological electronic stimulator (Jinan, Shan- 
dong, China). The day before the experiment, the mice with GTCS 
induced by MES were prepared for the tests of the anticonvulsant 
effects of the compounds. Mice were separately administered with 
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the testing compounds or vehicle in 5% dimethylsulfoxide/95% 
(0.2%) hydroxymethyl cellulose by intraperitoneal injection 30 
min before electroshock application in the mouse; the compound 
retigabine was used as the positive control. The incidence of hind 
limb extension was recorded, and if extension did not occur in the 
experiment, the animal was thought to be protected against the 
MES-induced seizure. For the PTZ-induced model, compounds 
were prepared in the same way as those in the MES-induced model 
and administered to male KM mice by intraperitoneal injection in 
a volume of 0.2 ml/ 10 g 30 min before administration of 100 mg/ 
kg PTZ sc. The latency to clonic seizure was recorded in seconds, 
and the incidence of GTCS and mortality rate caused by PTZ were 
both measured in this test. 

Test compounds 

Synthesis and characterization of YGPl (N-(6-Chloropyridin- 
3-yl)benzamide) To a solution of benzoic acid (0.36 mmol) in 2 
ml THE at 0 °C, A/,7V'-dicyclohexylcarbodiimide (DCC) was added 
with stirring for 30 min. 6-chloropyridin-3-amine (0.3 mmol) was 
then added to the mixture followed by addition of dimethyl ami- 
nopyridine (DMAP, 0.06 mmol) at room temperature with stirring 
overnight. After evaporating the volatiles, the mixture was parti- 
tioned by CH2CI2 (2 X 20 ml) and saturated with NaHCOj solution 
(20 ml), dried over anhydrous Na2S04 and concentrated in vacuo. 
The residue was purified by column chromatography to give the 
compound YGPl with 82% yield. HPLC (Agilent 1100; column/ 
Extend-C18, 5 fim, 4.6 x 150 mm; methanol:water/60:40; 0.6 ml/ 
min), > 98%; 'H NMR (400 MHz, CD3OD) 6 8.73-8.72 (m, IH), 
8.26-7.23 (m, IH), 7.96-7.93 (m, 2H), 7.62-7.58 (m, IH), 7.54-7.46 
(m, 2H), 7.44 (d, J= 0.8 Hz, IH); "C NMR (125 MHz, CD3OD) 
5 167.54, 145.11, 141.35, 135.28, 134.05, 131.95, 131.26, 128.34, 
127.35, 124.06; MS (ESI) m/z: 233 [M+H]*; HRMS (ESI) 
255.0301 CalcdforC12H9N20ClNa, found 255.031 1 [M+Na]*. 

Synthesis and characterization of YGPl (4-Fluoro-N-(pyridin- 
3-yl)benzaniide) To a solution of 4-fluorobenzoic acid (0.36 
mmol) in 2 ml THE at 0 °C, DCC was added with stirring for 30 
min. Then, pyridin-3 -amine (0.3 mmol) was added to the mixture 
followed by addition of DMAP (0.06 mmol) at room temperature 
with stirring overnight. After evaporating the volatiles, the mixture 
was partitioned by CHjCU (2 x 20 ml) and saturated with NaHCO, 
solution (20 ml), dried over anhydrous Na2S04 and concentrated 
in vacuo. The residue was purified by column chromatography to 
give the compound YGP2 with 85% yield. HPLC (Agilent 1100; 
column/Extend-C18, 5 \im, 4.6 x 150 mm; methanol:water/60:40; 
0.6 ml/min), > 99%; 'H NMR (CDCI3, 400 MHz): 8 8.65(d, IH), 
8.38(dd, IH), 8.28(m, IH), 8.10 (s, IH), 7.92 (m, 2H), 7.33(t, IH), 
7.16(dd, 2H); "C NMR (125 MHz, CDCI3+CD3OD): 5 166.3, 
164.1, 144.5, 141.4, 135.9, 130.5, 130.1 x 2, 128.6, 124.0, 115.6 x 
2. MS (ESI) m/z: 217 [M+H]"; HRMS (ESI) 239.0597 Calcd for 
C12H9N20FNa, found 239.0592 [M+Na]*. 

Structure characterization Compounds purchased from SPECS 
were purified before bioassay. The purity of the compounds listed 
in Figure 4B has been determined by HPLC to be >95%. 

Homology modeling of the transmembrane domain ofKCNQl 
KCNQ2 structures were modeled based on the structural infor- 
mation of the open-state Kvl.2 channel [29-31] by using Discov- 
ery studio 2.6 (Accelrys Inc., San Diego, CA, USA). Multiple se- 



quence alignment was generated by using CLUSTALW web server 
(http://www.ebi.ac.uk/Tools/msa/clustalw2) [48]. The highly con- 
served residues were used to guide the alignment (Supplementary 
information. Figure SIA), and the alignment was manually adjust- 
ed based on the transmembrane regions in the Kvl.2 structures [30, 
31, 49]. The result of the sequence alignment is shown in Supple- 
mentary infomation. Figure SIA. The S4-S5 linker, S5, and S6 
regions of Kvl.2 and KCNQ channels are highly conserved, which 
directly provides a reliable basis for mapping the positions of these 
regions of KCNQ channels based on the structure of Kvl.2 chan- 
nel. The secondary structure and transmembrane domain were pre- 
dicted by I-TASSER server [50] to guide the modeling of the Sl- 
S4 helix. Although residues are not so highly conserved in S1-S4 
segments as those in S5 and S6 segments, the conserved residues, 
EO, El, D3, Rl, R2, R3, R4 and R5, provided enough informa- 
tion to locate the positions of S 1 -S4 (Supplementary information, 
Figure SIA). The conserved phenylalanine of S2 also helped to 
map the position of S2. As a result, the positions of S1-S4 could 
be accurately mapped. The modeled structure of KCNQ2 channels 
underwent a 20 000-step energy minimization in AMBER 10 [51]. 

MD simulation 

MD simulations were performed on two systems — the 
KCNQ2 tetramer and the tetramer of the ztz240/KCNQ2 complex. 
MD simulation of KCNQ2 was started from the homology model 
of KCNQ2, and simulation of the complex started from the dock- 
ing model of ztz240 to the KCNQ2 homology model. The starting 
structures of the three systems were embedded in a patch of POPC 
bilayer generated by using VMD [52], and then solvated by TIP3P 
water molecules [53]. MD simulation was carried out by using 
GROMACS 4.53 [54]. The CHARMM27 force field with CMAP 
correction was applied [55-57]; ztz240 parameters were gener- 
ated by SWISSPARAM web server (http://www.swissparam.ch/) 
[58]. Na* and CI ions were added to maintain neutralization of the 
system. The NPT ensemble and 3D periodic boundary conditions 
were applied in the MD simulation, with the time step of 2 fs. The 
LINCS algorithm [59] was used to constrain the bond lengths. A 
constant pressure of 1 bar was applied with a coupling constant 
of 1.0 ps using Berendsen algorithm [60]. The v-rescale method 
[61] was employed to keep the temperature constant (300 K), with 
a coupling constant of 1.0 ps. The Particle-Mesh Ewald (PME) 
method [62] was used to handle the long-range electrostatics with 
a real-space cut-off of 1 .2 nm. The bilayer was first equilibrated to 
10 ns with the protein and ztz240 fixed. The systems of KCNQ2 
tetramer and ztz240/KCNQ2 complex were submitted to 50-ns and 
200-ns production runs, respectively. 

Docking and virtual screening 

Docking was perforaied by using Glide (Schrodinger, LLC, 
New York, NY, USA). Receptor was prepared by using the Protein 
Preparation and Grid Preparation tools in the Schrodinger Maestro 
interface. Hydrogen orientations were optimized in protein prepa- 
ration module. The docking box was set to cover the assumed 
binding pocket and centered at El 30 and R207, which were sup- 
posed to be the middle of the potential pocket. The default settings 
were adopted for the cut-off, neutralization, scaling and dimension 
of the binding pocket. According to the primary mutagenesis and 
electrophysiological results, we noticed that ztz240 may access 
F137. Homology modeling could not predict proper conformations 
for the residues interacting with a ligand. Therefore, we applied 
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the Induced Fit module encoded in the Schrodinger modeling 
package (Schrodinger, LLC). Inducted Fit is a flexible docking 
method that may consider the flexibilities of both the side chains 
of protein resides and small molecules during docking simulation. 
F137 and residues around F137 in the central part of VSD, E130, 
R207, S179, Q204, T133 and 1134 were treated as flexible resi- 
dues. 

The modules of High Throughput Virtual Screening (HTVS) 
and Standard Precision (SP) were used as filters sequentially in 
the virtual screening. From the commercially available SPECS 
database, ~20 000 compounds were screened and the top 700 com- 
pounds with highest G-score values were inspected. Compounds 
with the docking positions showing high-energy conformations, 
such as unreasonable angles and torsions caused by the steric ef- 
fect, were excluded. Among the remaining compounds, molecules 
with molecular weight > 500 and those with large subgroups out- 
side the pocket were also removed. The absorptions of the remain- 
ing potential hits were predicted by using Pipeline Pilot (Accelrys). 
AlogP < 5 and Absorption Level = 0 were taken as filters. Finally, 
25 compounds were selected as hits for bioassay. 

Volume calculation 

The sizes of the pockets in the KCNQ2 channel were calculated 
by using the program of POVME [63]. The sizes of ligands and 
the side chain of arginine were measured by the web server Molin- 
spiration (http://www.molinspiration.com/). 

Data and statistical analysis 

Patch-clamp data were processed using Clampfit 9.2 (Molecular 
Devices), and then analyzed in Graphpad Prism 4 (GraphPad Soft- 
ware, San Diego, CA, USA). Voltage-dependent activation curves 
were fit with the Boltzmann equation: G = G^,,, + (Gn,„ - G„i„)/ 
(1 + exp(K- Vy^lS)), where G„„ is the maximum conductance, 
G^i„ is the minimum conductance, K,,, is the voltage for reaching 
50% of maximum conductance and S is the slope factor. Dose- 
response curves were fit with the Hill equation: E = E„,J( 1 + (EC50/ 
C)""), where EC50 is the drug concentration producing half of the 
maximum response and P is the Hill coefficient. The deactivation 
traces were fit with exponential equations containing one or two 
components using Clampfit 9.2. Data are presented as means ± 
SEM. Significance was estimated using paired two-tailed Student's 
;-tests. 
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